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Nomenclature 


V S7<? 


A± 

= 

Rx ± - 2(o/6) 2 

a,b 

= 

plate length and width, see Fig. 1 

5± 

= 

(fi ± 7rU) + (a/bYR y± - ( a/by 

At 

= 

plate flexural stiffness 

h± 

= 

thickness of plate 

k 

= 

elastic spring constant 

l 

= 

lateral aerodynamic load 

M 

= 

Mach number 

Nx a 

N Vi = 

midplane force intensities, positive in compression 

7 

= 

dynamic pressure, p{/ 2 /2 

Rx± 

= 

a 2 N I± /ir 2 Z> ± 

Ry t 

= 

o ! iV y± / tt 2 D± 

S± 

= 

spring stiffness parameter, ka'/ir'D^. 

t 

= 

time 

u 

= 

freestream velocity 

w± 

= 

lateral deflection of plate 

x,y 

= 

Cartesian coordinates, see Fig. 1 

0 

= 

(Af 2 - l) 1 ' 2 

7± 

= 

mass density of plate 

X 

= 

2?a 3 //3D + 

P 

= 

mass density of air 

u> 

= 

circular frequency 

«± 2 

= 

ui 2 a i y±h±/D± 

+ , — 

= 

subscripts refer to upper and lower plate, respec- 


Introduction 

T HE flutter behavior of a structural configuration consist- 
ing of two rectangular, simply supported, parallel plates 
laterally connected by many closely spaced linear springs is 
investigated. The configuration analyzed is shown in Fig. 1. 
The upper plate has air flowing at supersonic speed over the 
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Fig. 1 Configuration and coordinate system. 


upper surface, and both plates are subjected to midplane 
loadings. 

This configuration is an idealization of a micrometeoroid 
bumper that is attached to a primary structure by a light, 
soft filler material. The aeroelastic behavior of such a con- 
figuration may be important in the design of. structural com- 
ponents of a manned space, station which are exposed to an 
airstream during launch. 


Analysis 

The equilibrium equations and appropriate boundary con- 
ditions are 

n , ,, C ) 2 u> + d*w + d‘ l w + 

D + VH» + + at ,+ + at v+ + y + h + — + 

k(w+ - w_) = l(x,y,t) (1) 
^ d 2 w>_ , , T d 2 u>_ , , dhu- . 

+ A *- -^r + - dlf +y- h - -w + 

k(w _ — w + ) = 0 (2) 

\ w ± (x,0 ,t) = w ± (x,b,t ) = w ± (0,y,t) = 

w ± (a,y,t) = Ol 


d*wj. 

dy 2 


(*,0,0 


^ (xM . 


dy 2 


d*w. 
dx 2 

dx 2 


(0,y,<) = 

(a,y,t) = 0| 


(3) 


R„_ = 0 
ok = 1 
R yi " 0 



Fig. 2 Flutter boundary for R x - — 0. 


r„. = R , 

o h = 1 



Fig. 3 Flutter boundary for R x + = R x - = R x . 


The frequency w is, in general, complex; however, attention 
is directed primarily to real values for which the motion is 
harmonic. 

When Eq. (4) is substituted into Eqs. (1) and (2) and the 
Galerkin procedure is used, the following equation is obtained : 


\l — A. 


B 4- -j- S-i 


S+S - 

1 - A_ - B_ + S- 


§A 

3 x 4 


8 

3 x 4 


16 - 4A + - B+ + S + - 


S+S- 

16 - 4 A_ - B_ + S_ 


= 0 


(5) 


where l(x,y,t) is the lateral load per unit area due to aerody- 
namic pressure. For static strip theory the lateral load is 
given bv the simple Ackeret value l{x,y,t) = —(2 q/ff) X 
(dw + /dx). 

A two-term Galerkin solution is pursued. Solutions that 
satisfy the boundary conditions for simply supported edges 
are assumed as follows: 

w±(x,y,t) = [c u± sin^™^ sin (Y) + 

C 2 , ± sin( 2 -^)sm(^)]e 4 “' (4) 


Flutter occurs with the coalescence of two natural frequen- 
cies as the dynamic pressure parameter X increases (see Ref. 
1). The procedure is to solve Eq. (5) for X and maximize the 
resulting expression with respect to the frequency to obtain a 
critical value of the dynamic pressure parameter X„. In 
many cases, more than one critical value exists corresponding 
to the coalescence of different pairs of modes, and it is neces- 
sary to seek the lowest critical value to define a flutter bound- 
ary. 

In order to illustrate the general flutter characteristics 
exhibited by this configuration, calculations were made for the 
simplified case for which h + = A_, D+ = /)., R v± = 0, a/b 
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R ». - 0 
o b = ! 

S = 20 



Fig. 4 Flutter boundary for Rz+ = 0. 


= 1, and S + = S_ = S. Flutter boundaries were derived from 
Eq. (5) for several combinations of streamwise midplane 
loads in the two plates. It must be remembered that these 
boundaries are subject to the same limitations that are de- 
scribed in Ref. 1 ; in particular, the range of validity of the 
boundaries is limited by the buckling characteristics of the 
plates. 

Discussion 

For a single flat isotropic plate, the flutter boundary as de- 
termined from a two-term Galerkin solution is a linear func- 
tion of the midplane load in the streamwise direction. Be- 
cause of the coupling of the motions of the two plates, how- 
ever, the configuration analyzed herein exhibits entirely dif- 
ferent boundaries. Peaks and valleys occur in the boundaries 


because the system can be tuned by means of the midplane 
loads. 

Figures 2-4 present flutter boundaries (A„ vs R*) for 
square plates having various combinations of midplane loads 
and a spring parameter of S = 20. At the present time 
realistic values of S are not clearly defined; the value chosen 
(>S = 20) might be typical of a configuration with a very soft 
filler material. The flutter boundary for a single flat plate 
with the same physical properties as either the upper or the 
lower plate considered in the present analysis is also shown in 
each of the figures (see Ref. 1). 

Figure 2 is a plot of the boundary for the configuration 
when there is no load in the lower plate. This boundary be- 
comes asymptotic to the single plate boundary for large nega- 
tive values of R„ i.e., large midplane tension, as do all of the 
boundaries considered herein. 

Figure 3 is a plot of the boundary when the midplane loads 
are the same in each plate. Here the tuning effect can be very 
significant since it is possible to have a zero flutter speed with 
a tensile load, and a peak exists which is much higher than the 
corresponding value for the single plate. 

Figure 4 is a plot of the boundary when there is no load in 
the upper plate. This case is perhaps the most realistic com- 
bination of loading if the configuration is considered to be a 
micrometeoroid protection device. This boundary has 
characteristics very similar to the boundary in Fig. 2 except 
that here the tuning effect is more prevalent. In particular, 
for streamwise tension, a condition that can be expected from 
bending loads on a space vehicle during launch, the elastically 
supported plate is more prone to flutter than the single plate 
alone. 

The results of this analysis indicate that if a configuration 
similar to this one is used for applications where supersonic 
airflows are encountered, a very careful flutter analysis is in 
order to insure that undesirable flutter characteristics are not 
present. 
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